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Abstract 

In this paper, we consider a point to point full duplex (FD) MIMO communication system. We 
assume that each node is equipped with an arbitrary number of antennas which can be used for 
transmission or reception. With FD radios, bidirectional information exchange between two nodes can 
be achieved at the same time. In this paper we design bidirectional link selection schemes by selecting 
a pair of transmit and receive antenna at both ends for communications in each direction to maximize 
the weighted sum rate or minimize the weighted sum symbol error rate (SER). The optimal selection 
schemes require exhaustive search, so they are highly complex. To tackle this problem, we propose 
a Serial-Max selection algorithm, which approaches the exhaustive search methods with much lower 
complexity. In the Serial-Max method, the antenna pairs with maximum “obtainable SINR” at both 
ends are selected in a two-step serial way. The performance of the proposed Serial-Max method is 
analyzed, and the closed-form expressions of the average weighted sum rate and the weighted sum 
SER are derived. The analysis is validated by simulations. Both analytical and simulation results show 
that as the number of antennas increases, the Serial-Max method approaches the performance of the 
exhaustive-search schemes in terms of sum rate and sum SER. 
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bidirectional link selection, full duplex, Serial-Max selection method. 

I. Introduction 

Current wireless communications systems typically exploit half duplex (HD) transmission. This 
is because for many years the full duplex (FD) transmission has been considered impractical. The 
signal leakage from the local output to input in the FD radio, referred to as the self interference, 
may overwhelm the receiver, thus making it impossible to extract the desired signals. Very 
recently, there has been a significant process in self-interference suppression in FD radios. The 
passive suppression methods |Q]-[[8l design the antennas using a combination of path loss, cross 
polarization and antenna directionality, while the active approaches 0], S, fl6], ifSl- lfTOl exploit 
the knowledge of self interference in cancelation in the analog or digital domain. The residual 
interference still exists and can be modeled as Rayleigh fading 0, llT6l . Ifl7l when the direct 
link is effectively suppressed. 

These suppression techniques can significantly reduce the self-interference, which has made 
FD radios practically feasible in the near future. This significant progress in FD has recently 
inspired some very interesting work on FD signal processing. In [18], ifTTl — lfT5l . the theoretical 
limits of point to point bidirectional FD have been investigated by taking into account the residual 
self-interference after suppression. In [0, the authors derived lower and upper bounds of the 
achievable sum rate of bidirectional FD communications, and proposed a transmission scheme to 
maximize the lower bound. In [fill , the achievable sum rate of bidirectional FD MIMO systems 
was analyzed and compared to the conventional HD MIMO systems over a spatial correlated 
channel. The ergodic capacity of bidirectional FD transmission using one transmit antenna and 
multiple receive antennas in the presence of channel estimation error has been derived in [jl2|. a 
FD antenna mode selection scheme was investigated in llT3l for a simple 2x2 MIMO system, 
where each antenna is either configured as the transmit or receive antenna mode. In llT4l . llT5ll . 
the suboptimal and optimal dynamic power allocation schemes were developed based on the 
sum rate maximization criterion. 

In this paper, we consider a general FD MIMO system with Na and N B antennas equipped 
at two nodes. Such an FD MIMO system will create NaN b possible links between the two 
FD MIMO nodes, with one possible link representing the channel from a transmit antenna of a 
node to a receive antenna of the other node. Since FD radios enable simultaneous bidirectional 
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information exchanges between two FD MIMO nodes, a fundamental question arisen in such 
a system is how to select the link for each direction to optimize the system performance. In 
this paper we consider two performance metrics, weighted sum rate maximization and weighted 
sum symbol error rate (SER) minimization. The optimal approach requires the exhaustive search 
from all possible antenna links, however, as the number of antennas increases, such a brute-force 
search bears very high complexity in selection process. 

To resolve this issue, in this paper we propose a simple Serial-Max selection algorithm by 
selecting the link with optimal performance for each direction in a two-step serial way, which 
can achieve asymptotically optimal performance. By using the law of total probability and order 
statics, the probability distribution functions of the two selected links are calculated, based on 
which, the closed-form expressions on average weighted sum rate and sum SER are derived. 
We show that the Serial-Max method approaches the brute-force search method in terms of the 
average weighted sum rate and sum SER as the number of antennas increases. The theoretical 
results are verified by Monte-Carlo simulations. 

The rest of this paper is organized as follows: Section II introduces the system model. The 
proposed Serial-Max selection algorithm is presented in Section III. Section IV analyzes the 
performance of the Serial-Max method, including the average weighted sum rate and sum SER. 
Simulation results are provided in Section V. In Section VI, we draw the main conclusions. 

II. System Model 

In this paper, we consider a bidirectional communication scenario between a pair of FD 
transceivers, node A and B, as illustrated in Fig. Q} where node A and B are equipped with 
Na and Nb antennas, respectively. Both nodes use the same frequency band at the same time 
for FD operation. Each node employs only one transmit and one receive RF chains, and any 
antenna can be configured to connect either the transmit or receive RF chain. In the proposed 
simultaneous bidirectional link selection (SBLS) scheme, two antenna links are selected for 
simultaneous bidirectional communication by selecting a pair of transmit and receive antennas 
at both ends. Within each antenna pair, one antenna is selected for transmission and one is for 
reception. 

’“Optimal” in this paper means that this scheme can achieve the optimal performance under practical constraint, that only 
the distribution of the residual interference rather than the instantaneous one can be obtained at each node. 
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Node A 


Node B 







Fig. 1. Full duplex MIMO systems with simultaneous bidirectional link selection 


We assume that the links between the two nodes are reciprocal and subject to indepen¬ 
dent Rayleigh fading, and together denoted by a Na x Nb channel matrix H = \h lJ ], i = 
1, 2, Na, j = 1, 2, Nb- The entry h l] represents the fading coefficient from the i-th antenna 
at node A to the j-th antenna at node B, and it follows the circularly symmetric Gaussian 
distribution CJV(0, a^). All the possible communication channels are assumed to follow the 
non-selective independent block fading, where the channel coefficients remain constant during 
a time slot, and vary from one to another independently. In the beginning of each time slot, all 
the possible communication links can be estimated perfectly. 

Since the two nodes are operated in the FD mode, there exists self interference caused by the 
signal leakage from the transmit antenna to the receive antenna at the same node. We assume the 
passive propagation suppression and active analog/digital cancelation techniques are employed to 
cancel the self interference. As indicated in [fl6ll . ifTTl . the direct link of the self interference can 
be effectively suppressed, and the residual self interference can be approximated to follow the 
Rayleigh distribution ll20ll . lf2Ti . In this paper, due to practical constraints, we assume that only 
the distribution of the residual interference, including the mean and variance, can be obtained 
and used in the selection process. Therefore, the selection process is based on the “obtainable 
SINR”. This is equivalent to that based on the instantaneous SNR, which will be later shown 
that the former is a scaled version of the latter. 

In the SBLS, the bidirectional links, i.e., (It- Jr), from the I T -\h antenna at node A to the 
J R -th antenna at node B, and (J T , I R ), from the J T -\h antenna at node B to the Ir -th antenna 
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at node A, are selected. The received signals at node A and B, denoted by y A and y B , can be 
expressed as 

Va = s/P t h jTlR x B + \[Pth^x A + n A (1) 

Vb = \[Pth lTjR x A + y/P t hfx B + n B , 

where h IrjR and h J ' r Ir denote the links corresponding to the selected antenna pair from node 
A to B and that from node B to A. P t is the transmit power at each node. The second term 
denotes the residual interference, at nodes A and B, respectively. We assume that both residual 
interference links, h A ! and h///, are subject to Rayleigh fading with zero mean and variance 
a BI . The AWGN at nodes A and B are denoted by n A and n B , which both follow CM {{), a 3 )- 
Without the knowledge of the instantaneous residual interference, the selection procedure is 
based on the “obtainable SINR” matrix T = [Y J ]. The entry in the “obtainable SINR” matrix 

• • i j . . p I u i j I 2 

is defined as Y ’ 3 = p-, where Ys = * a ' 3 ’ s the instantaneous SNR, and A* is the average 

p t < j2 

INR. Ys follows an exponential distribution with mean A., = — Y- and the average INR is given 

p.g. 2 . . 

by A i = r]X s = RI , where // denotes the cancelation ability. The “obtainable SINR” Y ' 3 is a 
scaled version of the instantaneous SNR Y'Y because we assume that all the links have the same 
average INR 

In the following analysis in Section IV, we first calculate the instantaneous performance based 
on the instantaneous residual interference, and then average it with the distribution of fading 

ip Jr 

channel and residual interference channel. We define the instantaneous SINR 7 AB = Yrrrr and 

1 Tb 

7 ba = liti . | for the two selected links, where the instantaneous INR 7^ 7 and yjj 7 are both 
exponential random variables with mean A,. 

III. Simultaneous Bidirectional Link Selection (SBLS) 

In this section, we first introduce two optimal SBLS approaches to maximize the weighted sum 
rate and minimize the weighted sum SER, respectively, based on the “obtainable SINR” matrix, 
or equivalently the SNR matrix. Then, a low-complexity method which achieves asymptotically 
optimal performance, referred to as Serial-Max method, is proposed. 

A. SBLS based on Weighted Sum Rate Maximization Criterion (Max-WSR) 

In this subsection, we describe the SBLS based on Weighted Sum Rate Maximization criterion 
(Max-WSR) under the Gaussian input assumption. In this criterion, two communication links 
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{(It, Jr), ( Ir, Jt )} from the It- th transmit antenna at node A to the Jr- th receive antenna at 
node B and the Jt- th transmit antenna at node B to the Ir -th receive antenna at node A, are 
selected to maximize the weighted bidirectional sum rate 

{(I T , Jr), (Ir, Jt)} = argmax {wR (Y tJr ) + (1 - w)R (f*) } , (2) 

1 <^r. A 

1 <jt ,jr<N B 

where R("f) = log 2 [1 + 7 ] denotes the rate under the “obtainable SINR” 7 = (A 7 ^ , and 0 < 
w < 1 is the given weight of the transmission from node A to B, depending on the rate 
requirement or quality of service (QoS) of each user. 

B. SBLS Based on Weighted Sum SER Minimization (Min-WSER) 

In the SBLS based on weighted sum SER minimization criterion (Min-WSER) under the 
assumption that the input signal is modulated with finite constellations, the bidirectional antenna 
links are selected to minimize the weighted sum SER 

{(I t ,Jr),(Ir,Jt)}= argmin {wSER (f rjt )+(l - w)SER } , (3) 

l<it,i r <N A 

}<jt,jr<N B 

where SER ( 7 ) = aQ (\JJy) represents the SER under the “obtainable SINR” 7 = • Q(') 

is the Gaussian Q-function ll22)l . and (a,/3) is a pair of constants determined by the modulation 
format, e.g., a = l,/3 = 2 for BPSK modulation. 

C. The Proposed Serial-Max SBLS 

The aforementioned SBLS schemes for the Max-WSR or Min-WSER criteria both require the 
brute-force search in order to find the optimal antenna pairs. This will become highly complex 
in selection process as the number of antennas increases. In this subsection, we introduce a 
low-complexity selection method, referred to as Serial-Max method, which selects the antenna 
pairs with maximum “obtainable SINR”, or equivalently, the maximum SNR, in a two-step serial 
way. 

In the first step, the best link with the maximum “obtainable SINR” is selected 

(Ji,Ji)= argmax {7 * m }. (4) 

1<H<JW<J1<JVB 
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We use 7 1 st = 7 /lJl to denote the “obtainable SINR” of the selected link (R. JR), and 7 to 
denote the SNR of this link. By removing the / r th column and J \-th row from the “obtainable 
SINR” matrix T (and the corresponding SNR matrix r s ), we can obtain a (N& — 1) x (N B — 1) 
submatrix T' (and a corresponding pruned SNR matrix r' s ). 

In the second step, the link (J 2 , J 2 ) with the maximum SNR is then selected from the pruned 
submatrix T' 

(J 2 ,J 2 )= argmax {Y rJt } ■ (5) 

1 < i 2 < Na , 1 <32 < Nb 
*2 7^1,72 7^1 

We use 7 2nd = Y 2j2 to denote the “obtainable SINR” of the second selected link by the Serial- 
Max method, which is the maximum element of T', and we use 7 2nd to denote the SNR of this 
link, which is the maximum element of the corresponding pruned SNR matrix T(. 

To maximize the weighted sum rate or minimize the weighted sum SER, the time-shared 
scheme can be employed, which allocates a fraction a of the time to use the best link for node 
yl’s transmission, and the rest 1 — a to use the best link for node B’s transmission. Given w, the 
weighted maximization sum rate problem can be solved by calculating the allocation fraction a. 
We have 

R = w[aR ( 7 ls4 ) + (1 - a) R (j 2nd )] + (1 - w) [aR (7 2nd ) + (1 - a) R ( 7 lst )] . (6) 

To maximize R, we first calculate the derivative 

^ (2»- - 1) [R ( 7 lst ) - R (' l 2nd )} • (7) 

It is obvious that R ( / y lst ) — R (^/ 2nd ) is positive, therefore the allocation factor a depends on 

wether w > 0.5 or not. Specifically, we have 

' 

a = 1, w > 0.5 

< a = 0, w < 0.5 • ( 8 ) 

0 < a < 1, w = 0.5 

Then, the weighted sum rate can be rewritten as 

R = max(w, 1 — w)R ( 7 lst ) + min(u>, 1 — w)R ( 7 2nd ), (9) 

The weighted sum SER minimization problem can be solved in a similar way, and we have 

SER = max(tu, 1 - w)SER ( 7 lst ) + min(u;, 1 - w)SER {^ 2nd ). (10) 
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It is shown from © and (fTOl) that the best link selected in the Serial-Max method will be 
allocated with the greater weight in the weighted sum rate and SER expression. 

Regarding the performance of the Serial-Max method, we have the following Lemma. 

Lemma 1: The Serial-Max method can achieve the optimal weighted sum rate and sum SER 
performances simultaneously, if the “obtainable SINR” of the second selected link r y 2nd is the 
second or third largest element of the “obtainable SINR” matrix T. 

Proof: We use R, and SERi to denote the respective rate and SER of the link with the /-1h 
largest “obtainable SINR”. If the second selected link corresponds to the second largest element 
in the “obtainable SINR” matrix T, the weighted sum rate is given by R = max(w, 1 — w)Ri + 
min(u’, 1 — ui'jRo. It is obvious that in this case the Serial-Max method can achieve the optimal 
performance in terms of weighted sum rate. Similarly, it can be easily proved that the Serial-Max 
method can achieve the optimal performance of SER. 

Recall that the pruned matrix T' is obtained by removing the row and column where the 
largest element, i.e., the first selected link, is located. Meanwhile, the second selected link is 
the largest element of the pruned matrix T'. Therefore, if the second selected link corresponds 
to the third largest element in the original “obtainable SINR” matrix, it implies that the second 
largest element in T is removed in the aforementioned manipulation. In other words, the two 
links associated with the first and second largest elements in T share the same antenna, which 
can not be selected simultaneously. In this case, the largest and third largest elements is the best 
option for the two selected links. Therefore, the Serial-Max method which selects the largest 
and third largest elements in T achieves the optimal performances of weighted sum rate and 
sum SER. ■ 

Then, we have the following proposition. 

Proposition 1: The probability that the Serial-Max method does not achieve the same perfor¬ 
mance as the optimal methods, denoted by P not , is upper bounded by 

p „ (N a + N b -2){N a + N b -S) 

(N A N B - 1) (N A N B - 2) ' K ’ 

Proof: According to Lemma 1, P not is upper bounded by the probability, denoted by P^, 
that the “obtainable SINR” of the second selected link R lnd is not the second nor third largest 
elements of the “obtainable SINR” matrix T. This implies that both the second and third largest 
elements are in the same row or column as the largest element, and they are both removed in 
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the process of obtaining the pruned “obtainable SINR” matrix T'. Note that the elements of T 
are independent and identically distributed. Due to symmetry, each element of E has the same 
probability to be the i-th largest element. Thus, we have 

p2,3 _ {Na + Nb — 2) (N a + N b — 3) 

(N a N b -1)(N a N b -2) ' { } 

2 3 

Combining the fact that P not is no more than P n ’ ot , Proposition 1 can be proved. ■ 

It is shown from (fTTI) that the upper bound of P not decreases quadratically as the number of 
antennas increases, which implies the probability that the Serial-Max method selects the same 
pairs as the optimal one will increase, and thus approaches the optimal performance in terms of 
weighted sum rate and sum SER asymptotically. 

In addition to the asymptotically-optimal performance, the complexity of the Serial-Max 
method is much simpler than the exhaustive search approach, as shown in Table |I] For the Serial- 
Max method, in the first step, the maximum “obtainable SINR” is selected from a Na x N b 
matrix, and N A N B comparisons are required. Similarly, for the second step (N A — 1 )(N B — 1) 
comparisons are needed, and the Serial method overall needs 2N A N B —N A —N B + 1 comparisons. 
By contrary, the optimal method requires exhaustive search in order to find the optimal antenna 
pairs, leading to NaNb ( Na ~ 1 ')( Nb ~ 1 ') comparisons. Therefore, the proposed Serial-Max algorithm 
can approach the optimal algorithm with significantly reduced complexity. 

TABLE i 

Complexity comparison 



Optimal selection approach 

Serial-Max approach 

Complexity 

N A N B (N A -1)(,N B -1) 

2 

2NaNb — Na — Nb + 1 


IV. Performance Analysis 

The optimal selection methods are both very difficult to analyze, and thus, in this section we 
analyze the performance of the Serial-Max algorithm. It will be shown later in simulations that 
the Serial-Max method can achieve near-optimal performance in terms of average weighted sum 
rate and SER. 
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A. Probability Distributions of Two Selected Links 


To analyze the performances of the Serial-Max method, the distributions of the real instanta¬ 
neous SINR corresponding to the two selected links 7 ab and 7 BA are required. For simplicity 
and without loss of generality, we consider the case of w > 0.5 for the following analysis, which 
can be easily extended to w < 0.5. 

According to the description of the Serial-Max method in ©, the SNR of the first selected link, 
7 I st is the largest order statistic among N A N B i.i.d. exponentially distributed random variables 
7 l J. The corresponding link (f, -J\) is used for transmission from node A to B. We have the 
following Lemma. 

Lemma 2: The CDF of 7 ab can be given by 

N A N B 

F 7 > AB (X)= 

k =0 


n a n b 

k 


- 1 / 


krjx + 1 


(13) 


Proof: The derivation is given in Appendix A. ■ 

The second selected link in the Serial-Max method (J 2 , Jf) is used for transmission from node 
B to A. According to ©, we can obtain 

Lemma 3: The CDF expression of the instantaneous link SINR of the second selected link, 
i.e., 7 ba, is given by 


-^7 BA (‘^) 


Na+N b —1 

E 


k =1 


N a N b l 

£ £ (—1 )"W 

l=N A NB~k m =0 


e 


N A N B -l + m 
A „ X 


(. N a N b — l + m) rjx + 


1 ’ 


(14) 


where Hk,i,m is expressed as 


( N A N B —k— 1 \ (N a N b \ ( l ^ 
\N A +N B -k-lJ \ l ) \m> 
dk,l,m ~ / N a Nb — 1 \ 

\N a +N b -2) 

Proof: The derivation is given in Appendix B. 


(15) 


B. Average weighted Sum Rate 

Based on the CDF expressions of the two selected links y AB and 7 BA , in this section, 
the average weighted sum rate of the two links is obtained. Firstly, the instantaneous rate is 
calculated for given realizations of the selected communication channel and the corresponding 
self-interference channel. Then the average rate is obtained by averaging the instantaneous rate 
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with respect to the distributions of the channels. Under the Gaussian input assumption, the 
average rate of the link with SINR 7 can be obtained as 


R = E[log 2 (l + 7)] 

OO 

1 f l-F^x) 
In 2 J 1 + x 


(16) 


dx, 


where FRx) is the CDF of 7. 

Proposition 2: The average weighted sum rate of the Serial-Max method, denoted by R s ~ Max 
is given by 

R s ~ Max = w r ab + (! - w)R B a, (IV) 

where Rab an d Rba are the average rates of the two selected links, respectively, which can be 
expressed as 


1 77 (-1)777) 

Rab — —- 2 _^ 


In 2 ' 1 — kri 

k =0 ' 


e^s Ei 


1 \ k 

— ) - e>sEi 
T]\, 


and 


N A +N B — 1 


= E W+^)' 


(18) 


(19) 


where 

N A N B -1 l 

* = E Ett 


1 \ m 


1 — (NaN b — l + m) 71 

l=N A N B —k m=0 \ A 


k=\ 


e r > x ‘> Ei 


1 

V^s 


n a n b —z+m / NaNb — l + m 

- e ^ Ei 


A, 


and 


N a N b , , m +l 
^ (_ -1 J Hk,N A N B ,m 


m= 1 


1 — 


e Es Ei ( —— I — e E Ei I — ) 


r/A, 


KJ 


( 20 ) 

( 21 ) 


In addition, Ei(-) denotes the exponential integral function 

Proof: The derivation is given in Appendix C. ■ 

It is shown from (fl7l)-(l2TI) that the average weighted sum rate of the Serial-Max method 
R s ~ -Max j s on ]y determined by the average SNR A, s and cancelation ability //. When A s goes to 
infinity, we have Ei(e) —7 — ln(e) and e E ~ 1 + e, for e = 1/A S , where 7 ~ 0.5772 is Euler’s 
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constant [f28l . After some manipulations, R s Max converges to a rate ceiling 

fbS—Max \ ^ ( 1) (, k ) (fc 7 ?) , 1 ~ w \ ^ \ ^ (~1) dk,N A N B ,m hi (™?) 

A s —>oo j ri 2 2.^ 1 — kri In 2 ^ ^ 1 — rmi 

k =0 ' fc=l m=l 

( 22 ) 


1 — w 
In 2 


A^+Afg — 1 NaNb~ 1 l 

E E E 

fc=l l=N A N B -km=0 


(-1 ) m Hk,i,m In [{N a N b -l + m)rj] 
1 - ( N a N b - l + m) r) 


C. Average weighted sum SER 


In this section, we analyze the average weighted sum SER of the Serial-Max method. For the 
SER analysis, we assume that the input signal is modulated with a finite constellation. Though 
the finite constellations like BPSK are used, the current self-interference cancellation technique 
can reduce the self-interference near to the noise level, as shown in ifTTTl . f[24l . It has thus been 
commonly assumed in many existing papers that the residual self-interference after cancellation 
follow the Rayleigh distributions lfT 6 l . ||20|| . We also adopt this assumption in this paper. Firstly, 
the SER is calculated for a given set of channel realizations, similar to the weighted sum rate 
analysis. Then, the SER is averaged over the communication and self-interference channels. The 
average SER of the link with SINR 7 can be written as 


SER = oE 


q(v% : 


(23) 


Q x 

Oiy/]3 f F^(x)e~~^ 

2 ^ 2 n J 


dx. 


Q(-) is the Gaussian Q-Function |{22| . 

In addition, if the first-order expansion of the PDF of 7 is expressed as 


(x 


N 


fj( x ) = TETT + °( xN+£ )i 


X N + 


the asymptotic SER can be obtained as 


SER = 


2 N a£T (N + |) 


0F(7V+1)(/3A ) jv+1 + ° U (7V+1) / ' 


1 \ 


Proposition 3: The average weighted sum SER of the Serial-Max method, denoted by SER 


(24) 


(25) 


S—Max 


IS 


SER 


•S—Max 


= wSERab + (1 — w)SER B Ai 


(26) 
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where SER AB and SER BA are the average SER of the two selected links, and can be calculated 
as 

_ i (-i) k ( NA . NB ) ( R T\ i , fl 

(27) 


i nrz n a n b ( -1 \k(N A N B \ 

SER AB = i + £ ~ i-* - Q 


V2 


k =1 


y/krj 


2 P \ _^ + jl 
-| -I e»;w , 


rj\ s krj 


and 


where 


’ Na+Nb~ 1 


/To— 

ser ba = ^ Y (G 1 + G 2 ), 


V2 


(28) 


fc=i 


P 


N A N B — 1 l 

G = V V ^ 

,-fci^o ^ ' WV fl -Z+771)77 


1 \ m 

't) !~^k,l,r. 


l=N A N B - 


| g»7W ^2(iV A JV s -i+m)j7 ^ 

(29) 


and 


G 2 - Y 


N a Nb , s m 

(-i) Hh,N A N B , 


P 


1 _l_£_ 


-<3 t/—r- H- e* 

?7A s m?7 


(30) 


m=l v 

Proof: The derivation is given in Appendix D. ■ 

The SER performance converges to an error floor, when the average SNR \ s increases to 
infinity 


m — n aNb ( 1 A+l (N A N B \ 

SERZZr Al - z - { k ] 


^S—Mcix 


(1 — w) a^/pTT 

71 


y/kp 


Q| ,/A) e i 


(31) 


Na+Nb — 1 / \ m+l 

l -1 ) Vk,N A N B ,m Q 


k =1 


y/mrj 


P \ JL. 

I g 2mrj 

mrj 


/ -1 \ rn — N a +Nb~ 1 N a Nb l / -.\m+l 

(1 — Ml) a\/p7r y-v (~1) Pk,l,m 


V2 


x Q 


k =1 l=N A N B —k m=0 7 (N A N B l + m) r) 

p 

g2(N A N b ~l+m)ri 


P 


(N a N b -l + m)r] 

On the other hand, when 77 = 0, i.e. the self interference is perfectly canceled, we can further 
calculate the asymptotic SER of the Serial-Max method at high SNR. Firstly, the CDF expression 
of 7 ab can be rewritten as 


F iab( x ) = U-e * 


N A N B 


Then, the first order expansion of its corresponding PDF is given by 




x 


n a n b —1 


B~ 


A 


N A N B 


+ o [x 


NaNB~ 


)■ 


(32) 


(33) 
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Using (l25l) . the asymptotic SER of 7 AB can be obtained 


SER ab = 


u 1 


A 


N A N B 


+ o 


A 


N A N B 


(34) 


where 


u 1 = 


2 N A N B-i aT ( N a N b + |) 


/3 N aN B yTf 

Similarly, when the self interference cancelation is perfect, the CDF of 7 BA is rewritten as 

N A +N B — 1 NaN b 


(35) 


^\ba ( 3 ) 


t: yi 

fc=l l=N A N B -k 


_* {N A N B -l)x 

1 — e e , 


where 


A 1 iti — 


( N A N B -k-l \ (N a N b \ 
\N a +N B — k—l) \ l ) 

( n a n b —i \ 

\Na+N b -2J 


(36) 


(37) 


Then, the first order expansion of its corresponding PDF is calculated as 


fl BA ( X ) — N a +N b -1,(Na-1)(N b -1)(Na 1)(^Vb 1) 


X 


(N a -1)(N b -1)-1 


(Na-1)(Nb-1) 


A 


+ o (a: (ArA " 1)(7Vs_1)+e ) 

(38) 


Combing (|25T) . the asymptotic SER of 7 B a is 


where 


SERba ^n a -i)(n b - r) + 0 ^(JV A -1)(JV B -1) ) ’ 

2 N A N B - NA - NBaT ( {Na _ 1){Nb _ 1} + i) ( (JVa 5 A ^_ 1} ) 


(39) 

(40) 


Combining (l34l) and (l39l) . the asymptotic weighted sum SER with perfect interference cance¬ 
lation is expressed as 


SER S 


—Max 


(1 -W)U 2 f 1 \ 

^(tV A -l)(AT fl -l) y^{N A -l)(N B -l) J ■ 


(41) 


It is implied by (l4ll) that given 0 < w < 1 the diversity order of the Serial-Max method is 
(N a — 1)(Nb — 1) with 77 = 0. This is coincident with a simple deduction of the existing 
result f271 : The diversity order is determined by the worse link, i.e., the second best link which 
is selected from T 1 ' consisting of (N A — 1 )(N B — 1) i.i.d. elements. Moreover, the transmission 
direction with the greater weight will achieve the full diversity order of N A N B , and the other 
direction will achieve the diversity order of (N A — 1 )(N b — 1) for w ^ 0.5. On the condition 
that w = 0.5, a can be a arbitrary fraction. Then the best link can be arbitrarily allocated to 
each direction. Therefore, the diversity orders of the two directions are both (N A — l)(iV B — 1), 
and the achievable diversity is obviously ( N A — 1 )(N B — 1) for the weighted sum SER. 
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V. Simulation Results 

In this section, we provide the simulation results for our proposed SBLS methods to validate 
the previous analysis. For simplicity, we consider that Na = N B = N in the following 
simulations. 

A. Average Weighted Sum Rate 

Fig. [2] depicts the average weighted sum rate of the Max-WSR and Serial-Max methods versus 
SNR with different levels of self interference rj = 0.02, 0.05, 0.1 for N = 3 and w = 0.7. It can be 
seen that the weighted sum rate expression in (fTTT) perfectly matches with the simulation results. 
In addition, the weighted sum rate performance is limited by rate ceilings, which coincides 
with the preceding analysis in (l22l) . From the figure, we find that at low A. s , the weighted sum 
rate performance for different 77 is quite similar, because the weighted sum rate performance 
at low A s is SNR-limited. However, at large A s , the residual self interference will dominate 
the performance, and the performance is limited by the rate ceiling caused by the residual self 
interference. The figure also reveals that the Serial-Max method achieves almost the same average 
weighted sum rate as the Max-WSR method across all SNR regions. 

In Fig. El we illustrate the average weighted sum rate of the Max-WSR and Serial-Max 
methods for different numbers of antennas N = 3,4,5 where the self interference cancelation 
coefficient is 77 = 0.02. It can be observed that, for different numbers of antennas, the Serial-Max 
algorithm achieves almost the same average weighted sum rate as the Max-WSR one across all 
SNR region. We can also find from this figure that the average weighted sum rate increases with 
the number of antennas. 

B. Average Weighted Sum SER 

The following simulations of weighted sum SER performance are conducted with BPSK 
modulation. In Fig. the weighted sum SER performance of the Serial-Max method is provided 
for different 77 = 0,0.05,0.1,0.5, where large 77 means severe self interference whereas small 
77 means slight self interference level. Especially, 77 = 0 means perfect interference cancelation. 
This figure verifies the weighted sum SER expression given by Proposition 3. Based on the 
figure, it can be observed that the simulated SER performance for 77 = 0 tightly matches with the 
asymptotic one given by (l4TT) at high SNR, while the SER performances for 77 > 0 are constrained 
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Fig. 2. Weighted sum rate performances of the Max-WSR and Serial-Max methods, where w = 0.7, N = 3 and rj = 
0.02,0.05,0.1. 



Fig. 3. Weighted sum rate performances of the Max-WSR and Serial-Max methods, where w = 0.7, rj = 0.02 and IV = 3,4, 5. 
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Fig. 4. Weighted sum SER performances of the Serial-Max method, where N = 3 and r} = 0, 0.05, 0.1, 0.5. 

by error floors evaluated by OTI) at high SNR. It can also be seen that when the self interference 
is perfectly canceled, the diversity order of the Serial-Max method is (IV 4 — l)(N B — 1). However 
when the residual self interference exists, the Serial-Max method has a zero diversity order. 

Fig. [5] shows the average weighted sum SER for different numbers of antennas TV = 3,4, 5 
where the self interference level 77 = 0.1 is assumed. It shows that both the Min-WSER and 
Serial-Max methods are limited by error floors at high SNR due to the residual self interference. 
As the number of antennas increase, the SER performances of both methods including the error 
floor are improved. Moreover, the Serial-Max method performs closer to the Min-WSER scheme 
as N increases. 

Fig. [ 6 ] compares the simulated weighted sum SER performance of the Min-WSER and Serial- 
Max methods with different number of antennas. Combinations of different SNR and self 
interference levels (A s = 10, 15dB, 77 = 0.1, 0.2) are provided. It can be observed that the gaps 
between the Min-WSER and Serial-Max methods are reduced as N increases. It also shows 
that the SER performance of the Serial-Max method approaches the Min-WSER method when 
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Fig. 5. Weighted sum SER performances of the Min-WSER and Serial-Max methods, where 77 = 0.05 and N = 3,4, 5. 

the self interference is large or SNR is small. This is because in these cases these two factors 
dominate the SER performances of both methods. 

C. Computational complexity 

To compare the computational complexity of the optimal and Serial-Max methods, the number 
of required floating-point operations (flops) are presented in Fig. [7] It is clear that the optimal 
method has a high complexity due to the brute-force search operation, while the Serial-Max 
method can provide significant a complexity reduction, especially when the number of antennas 
is large. 


VI. Conclusions 

In this paper, we have proposed an opportunistic bidirectional link selection approach in FD 
MIMO systems. The optimal approach based on the “obtainable SINR", defined as the ratio 
of the instantaneous SNR and average INR plus one, requires exhaustive search, bearing high 
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Fig. 6. Simulated weighted sum SER performances of the Min-WSER and Serial-Max methods, where r) = 0.1, 0.2 and 
Pt = 10,15dB. 


complexity. To facilitate the selection process and reduce the computational complexity, a simple 
Serial-Max method with near optimal performances was proposed. The performance analysis was 
provide for the Serial-Max method, including the weighted sum rate and SER performances. It 
was shown that the proposed Serial-Max method approaches the respective weighted sum rate 
and weighted sum SER performances of the exhaustive search methods when the number of 
antennas increases. 


Appendix A: Proof of Lemma 2 

According to the description of the Serial-Max method, the CDF of 7 ab is expressed as 


F 1ab( X ) = Pr 


ryhJl 
7 B 1 + 1 


< X 


= Pr <x( 7 | / + l)), 


(42) 
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Fig. 7. Complexity comparison between the optimal method and the Serial-Max method. 


where 7 f lJl is the largest one of JV 4 x A r /? i.i.d. exponential-distributed random variables, and 
its CDF is given by ll23l 



— x 
e x » 


n a n b 


(43) 


7 IP follows exponential distribution with average A*. Substituting the CDF expression of 7 J 1 J| 
and 7 ^, (1421) can be calculated as 




0 L - 


1 — e 


°(rg' + i) 




-d 7 | 7 - 


Through some manipulations, we can obtain the CDF of 7 ab- 


(44) 


Appendix B: Proof of Lemma 3 


As aforementioned, we have 7 B a = 7s / (7a + 1)’ an d its CDF expression is 


F 1ba( x ) = Pr 


7. 


I2J2 


7 ^ + 1 


< x 


(45) 


= Pr (7f 2J2 < a; (7 a + l)) • 
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''/I' 2 ' 1 ' 2 is the maximum order statistic in the pruned matrix T' s consisting of (N A — 1) x (N B — 1) 
random variables following the exponential distribution. Since T' s is obtained by removing the 
maximum element and other N A + N B — 2 elements from T s , 7 J 2 J ' 2 could possibly be one of 
any (N A N B — k)- th order statistic of 1%, k = {1, 2..., N A + N B — 1}. As a result, the CDF of 
7^ 2j2 is given by 

N A +N B — 1 

F y y a {x)= ^ PkFi NANB ~ k \x), (46) 

k =1 

where is the probability that 7 721/2 is the (N A N B -k)- th order statistic in T s , and F^ NANB ~ k \x) 
represents the CDF of the (N A N B — /c)-th order statistic among N A N B variables. 

For n i.i.d. variables, each with CDF F(x), the CDF of the r-th largest order statistic can be 
written as 

F {r \x) = J2 1 - F(x)} n -\ (47) 

n—r ' ' 


Here, they are all exponential variables with average X s . Then, the CDF of the (N A N B — /c)-th 
order statistic can be expressed as 


n a n b 

F W A N B -k) {x)= £ 

l=N A N B -k 
n a n b l 

= E E 

l=N A N B —k m =0 


n a n b 

l 


1 — e t 


(48) 


N a N b 

l 


N a N D — Z+m 

/ a \m - A - x 

(—1) e 1 


On the other hand, if 7 721/2 is the (N A N B — /c)-1h order statistic in T s , it means that among the 
(N a + N b — 1) elements removed from the matrix F. s , there are k elements larger than / 'jl 2 ' J ‘ 2 . To 
satisfy that the largest element of T' s is the (/;: + 1) largest element in T s , we can first remove 
the first k largest elements in the matrix T s , and then, keeping the (/;: + 1) -1h largest element 
unremoved, we remove the other (N A + N B — 1 — k) elements randomly. Thus, 77 . is calculated 
as 

V n a +n b - 


Pk = 


( N A N B -k-l \ 
\N A +N B -k-lJ 


( N a N b — 1 \ ’ 
\N a +N b -2) 


(49) 


where the numerator implies that when we remove (N A + N B — 1 — k) elements randomly from 
(N a N b — k — 1) elements, there are possibilities, while the denominator means 

that when we remove (N A + N B — 2) elements (besides 7 f lJl , the largest one of T s ) randomly 
from ( N a N b — 1) elements, there are ( a a ^ a a E' 2 ) possibilities. 
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Finally, substituting (l48l) and (l49l) into (1461) . the CDF expression of y ^ 2 ’ 72 can be obtained as 

Nj,i+Ng — 1 N A N B l 

^ ^ N a N id —Z+m 

]T (-l) m ^ m e -x-*, (50) 

k =1 l=N A N B —km=0 

where ftk,i, m is expressed as 


ft k ,1,171 


( N A N B —k—l \/ N A N B \/l \ 
\N A +N B -k-l)\ l A m) 


N a N b —1 
Na+N b 


%) 


Then, substituting (1501) into (1451) and combining the CDF of 7 ^ we have 


°?N a +N b -1 n a n b l 

F i ba( x )= 

n k =1 l=N A N B -km=0 


- NANB x ~ l+rn X {lf+l) e 


RI 

.1A 




-dj 


it/ 


Through some manipulations, Lemma 3 can be obtained. 


(51) 


(52) 


Appendix C: Proof of Proposition 2 

Combining the CDF expression of 7 ab in (fOl) and (fl6l) . we can obtain 

n a n b 




— X 

e 


n a n b 


In 2 


EM) 


k =1 


fc+i (NaNb 
k 


( kr/x +!)(! + x) 


-R x 

e Xs 


dx 


1 — kr) 


kr]e~^ x 


(1 + x) krjx + 1 


dx. 


Using the integral Ei (x) = e x / 0 °° f—.dd Il22l . (1531) can be given by 


D 1 ™(-l) fe (™) 

Jd.AB = 1 —- 2 _^ 


In 2 


k =1 


1 — kr] 


1 / 1 \ _fc ( k 

e r ,\ 3 E x —— — e x ° Ex — 

vKJ \K 


On the other hand, the CDF of 7 B a in (fl4l) can be rewritten as 

n a +n b —1 

F 1ba( X ) =1 ~ (/1 + /2), 

k=\ 


where 


and 


n a n b —1 1 

h = E E Mr + A*, 

l=N A N B —k m =0 


N A N b —Z+m 

e J 


( NaN b — l + m) r]x + 1 ’ 


n a n b 

$2 = ^ (—l) m+ Vfc,/V A IV s ,m 
m=l 


— 7 —x 

e Xs 


mr]x + 1 


(53) 


(54) 


(55) 


(56) 


(57) 
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Then, substituting (l 55 l)-(l 57 T) into ([T6l) . we have 


Rba — 


1 


Na+N b — 1 


In 2 


E 

k=1 


/l + /2 

1 + X 


dx 


(58) 


After some manipulations, the average rate Rba in (fT9b can be obtained. Finally, combining the 
expression of Rab and Rba, Proposition 2 can be proved. 


Appendix D: Proof of Proposition 3 
Firstly, substituting the CDF of 7 ab in (fl 3 l) into ( 1231 ) . we have 


SERab = -~ 


1 av^ 


2 2a/2e 


NaNb /at at \ r 

e <-e +i CT)./(^ 


-($+£)* 


(krjx + 1) yfx 


dx. 


Applying Integration by substitution, the integral part is rewritten as 

f e. U + aJ x t = Ci f 2e U + aJ* 

dx — —— -—at. 


0 

2 2 


(krjx + 1) J ( krjt 2 + 1) t 

0 


Using the integral f 0 °° ^^rdx = Q {y/ 2 ab) f e a2ft2 


can be calculated as 


e (2 a) ) x 27 t 

=aa; = _ Q 


+ 1) a/A?// \ V krj 


2 (3 \ _^ + jl 

_|-I g ?7 As ' 2fe»j _ 


(59) 


( 60 ) 


( 61 ) 


0 \ / 

Substituting (l6ll) into (l59l) . and after some manipulation, (l27l) can be obtained. 

Then, the average SER of 7 ba in ( 1281 ) can be obtained by substituting (fl4l) into ([ 23 b and 
through the similar manipulations. Finally, combining the two average SER expressions, Propo¬ 
sition 3 can be proved. 
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